SABRes: in silico detection of drug resistance conferring mutations in subpopulations of SARS-CoV-2 genomes.

The emergence of resistance to antiviral drugs increasingly used to treat SARS-CoV-2 infections has been recognised as a significant threat to COVID-19 control. In addition, some SARS-CoV-2 variants of concern appear to be intrinsically resistant to several classes of these antiviral agents. Therefore, there is a critical need for rapid recognition of clinically relevant polymorphisms in SARS-CoV-2 genomes associated with significant reduction of drug activity in virus neutralisation experiments. Here we present SABRes, a bioinformatic tool, which leverages on expanding public datasets of SARS-CoV-2 genomes and allows detection of drug resistance mutations in consensus genomes as well as in viral subpopulations. We have applied SABRes to detect resistance-conferring mutations in 25,197 genomes generated over the course of the SARS-CoV-2 pandemic in Australia and identified 299 genomes containing resistance conferring mutations to the five antiviral therapeutics that retain effectiveness against currently circulating strains of SARS-CoV-2 - Sotrovimab, Bebtelovimab, Remdesivir, Nirmatrelvir and Molnupiravir. These genomes accounted for a 1.18% prevalence of resistant isolates discovered by SABRes, including 80 genomes with resistance conferring mutations found in viral subpopulations. Timely recognition of these mutations within subpopulations is critical as these mutations can provide an advantage under selective pressure and presents an important step forward in our ability to monitor SARS-CoV-2 drug resistance.

Guiding the design of SARS-CoV-2 genomic surveillance by estimating the resolution of outbreak detection

Genomic surveillance of SARS-CoV-2 has been essential to inform public health response to outbreaks. The high incidence of infection has resulted in a smaller proportion of cases undergoing whole genome sequencing due to finite resources. We present a framework for estimating the impact of reduced depths of genomic surveillance on the resolution of outbreaks, based on a clustering approach using pairwise genetic and temporal distances. We apply the framework to simulated outbreak data to show that outbreaks are detected less frequently when fewer cases are subjected to whole genome sequencing. The impact of sequencing fewer cases depends on the size of the outbreaks, and on the genetic and temporal similarity of the index cases of the outbreaks. We also apply the framework to an outbreak of the SARS-CoV-2 Delta variant in New South Wales, Australia. We find that the detection of clusters in the outbreak would have been delayed if fewer cases had been sequenced. Existing recommendations for genomic surveillance estimate the minimum number of cases to sequence in order to detect and monitor new virus variants, assuming representative sampling of cases. Our method instead measures the resolution of clustering, which is important for genomic epidemiology, and accommodates sampling biases.

Improved Neutralisation of the SARS-CoV-2 Omicron Variant following a Booster Dose of Pfizer-BioNTech (BNT162b2) COVID-19 Vaccine.

In late November 2021, the World Health Organization declared the SARS-CoV-2 lineage B.1.1.529 the fifth variant of concern, Omicron. This variant has acquired over 30 mutations in the spike protein (with 15 in the receptor-binding domain), raising concerns that Omicron could evade naturally acquired and vaccine-derived immunity. We utilized an authentic virus, multicycle neutralisation assay to demonstrate that sera collected one, three, and six months post-two doses of Pfizer-BioNTech BNT162b2 had a limited ability to neutralise SARS-CoV-2. However, four weeks after a third dose, neutralising antibody titres were boosted. Despite this increase, neutralising antibody titres were reduced fourfold for Omicron compared to lineage A.2.2 SARS-CoV-2.

Development and Validation of an In-House Real-Time Reverse-Transcriptase Polymerase Chain Reaction Assay for SARS-CoV-2 Omicron Lineage Subtyping between BA.1 and BA.2.

In order to rapidly differentiate sublineages BA.1 and BA.2 of the SARS-CoV-2 variant of concern Omicron, we developed a real-time reverse-transcriptase polymerase chain reaction to target the discriminatory spike protein deletion at amino acid position 69-70 (S:del69-70). Compared to the gold standard of whole genome sequencing, the candidate assay was 100% sensitive and 99.4% specific. Sublineage typing by RT-PCR can provide a rapid, high throughput and cost-effective method to enhance surveillance as well as potentially guiding treatment and infection control decisions.

SARS-CoV-2 Within-Host and in vitro Genomic Variability and Sub-Genomic RNA Levels Indicate Differences in Viral Expression Between Clinical Cohorts and in vitro Culture.

Background: Low frequency intrahost single nucleotide variants (iSNVs) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have been increasingly recognised as predictive indicators of positive selection. Particularly as growing numbers of SARS-CoV-2 variants of interest (VOI) and concern (VOC) emerge. However, the dynamics of subgenomic RNA (sgRNA) expression and its impact on genomic diversity and infection outcome remain poorly understood. This study aims to investigate and quantify iSNVs and sgRNA expression in single and longitudinally sampled cohorts over the course of mild and severe SARS-CoV-2 infection, benchmarked against an in vitro infection model. Methods: Two clinical cohorts of SARS-CoV-2 positive cases in New South Wales, Australia collected between March 2020 and August 2021 were sequenced. Longitudinal samples from cases hospitalised due to SARS-CoV-2 infection (severe) (n = 16) were analysed and compared with cases that presented with SARS-CoV-2 symptoms but were not hospitalised (mild) (n = 23). SARS-CoV-2 genomic diversity profiles were also examined from daily sampling of culture experiments for three SARS-CoV-2 variants (Lineage A, B.1.351, and B.1.617.2) cultured in VeroE6 C1008 cells (n = 33). Results: Intrahost single nucleotide variants were detected in 83% (19/23) of the mild cohort cases and 100% (16/16) of the severe cohort cases. SNP profiles remained relatively fixed over time, with an average of 1.66 SNPs gained or lost, and an average of 4.2 and 5.9 low frequency variants per patient were detected in severe and mild infection, respectively. sgRNA was detected in 100% (25/25) of the mild genomes and 92% (24/26) of the severe genomes. Total sgRNA expressed across all genes in the mild cohort was significantly higher than that of the severe cohort. Significantly higher expression levels were detected in the spike and the nucleocapsid genes. There was significantly less sgRNA detected in the culture dilutions than the clinical cohorts. Discussion and Conclusion: The positions and frequencies of iSNVs in the severe and mild infection cohorts were dynamic overtime, highlighting the importance of continual monitoring, particularly during community outbreaks where multiple SARS-CoV-2 variants may co-circulate. sgRNA levels can vary across patients and the overall level of sgRNA reads compared to genomic RNA can be less than 1%. The relative contribution of sgRNA to the severity of illness warrants further investigation given the level of variation between genomes. Further monitoring of sgRNAs will improve the understanding of SARS-CoV-2 evolution and the effectiveness of therapeutic and public health containment measures during the pandemic.

Co-infection with SARS-CoV-2 Omicron and Delta variants revealed by genomic surveillance.

Co-infections with different variants of SARS-CoV-2 are a key precursor to recombination events that are likely to drive SARS-CoV-2 evolution. Rapid identification of such co-infections is required to determine their frequency in the community, particularly in populations at-risk of severe COVID-19, which have already been identified as incubators for punctuated evolutionary events. However, limited data and tools are currently available to detect and characterise the SARS-CoV-2 co-infections associated with recognised variants of concern. Here we describe co-infection with the SARS-CoV-2 variants of concern Omicron and Delta in two epidemiologically unrelated adult patients with chronic kidney disease requiring maintenance haemodialysis. Both variants were co-circulating in the community at the time of detection. Genomic surveillance based on amplicon- and probe-based sequencing using short- and long-read technologies identified and quantified subpopulations of Delta and Omicron viruses in respiratory samples. These findings highlight the importance of integrated genomic surveillance in vulnerable populations and provide diagnostic pathways to recognise SARS-CoV-2 co-infection using genomic data.

Assessment of Inter-Laboratory Differences in SARS-CoV-2 Consensus Genome Assemblies between Public Health Laboratories in Australia.

Whole-genome sequencing of viral isolates is critical for informing transmission patterns and for the ongoing evolution of pathogens, especially during a pandemic. However, when genomes have low variability in the early stages of a pandemic, the impact of technical and/or sequencing errors increases. We quantitatively assessed inter-laboratory differences in consensus genome assemblies of 72 matched SARS-CoV-2-positive specimens sequenced at different laboratories in Sydney, Australia. Raw sequence data were assembled using two different bioinformatics pipelines in parallel, and resulting consensus genomes were compared to detect laboratory-specific differences. Matched genome sequences were predominantly concordant, with a median pairwise identity of 99.997%. Identified differences were predominantly driven by ambiguous site content. Ignoring these produced differences in only 2.3% (5/216) of pairwise comparisons, each differing by a single nucleotide. Matched samples were assigned the same Pango lineage in 98.2% (212/216) of pairwise comparisons, and were mostly assigned to the same phylogenetic clade. However, epidemiological inference based only on single nucleotide variant distances may lead to significant differences in the number of defined clusters if variant allele frequency thresholds for consensus genome generation differ between laboratories. These results underscore the need for a unified, best-practices approach to bioinformatics between laboratories working on a common outbreak problem.

Documenting elimination of co-circulating COVID-19 clusters using genomics in New South Wales, Australia.

OBJECTIVE: To adapt 'fishplots' to describe real-time evolution of SARS-CoV-2 genomic clusters. RESULTS: This novel analysis adapted the fishplot to depict the size and duration of circulating genomic clusters over time in New South Wales, Australia. It illuminated the effectiveness of interventions on the emergence, spread and eventual elimination of clusters and distilled genomic data into clear information to inform public health action.

Bondi and beyond. Lessons from three waves of COVID-19 from 2020.

OBJECTIVES: To describe local operational aspects of the coronavirus disease 2019 (COVID-19) response during the first three waves of outbreaks in New South Wales (NSW), Australia, which began in January, July and December 2020. Type of program or service: Public health outbreak response. METHODS: Narrative with epidemiological linking and genomic testing. RESULTS: Epidemiological linking and genomic testing found that during the first wave of COVID-19 in NSW, a large number of community transmissions went undetected because of limited testing for the virus and limited contact tracing of cases. The second wave of COVID-19 in NSW emerged following reintroduction from the second wave in Victoria, Australia in July 2020, and the third wave followed undetected introduction from overseas. By the second and third waves, cases could be more effectively detected and isolated through an increased ability to test and contact trace, and to rapidly genomic sequence severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) isolates, allowing most cases to be identified and epidemiologically linked. This greater certainty in understanding chains of transmission resulted in control of the outbreaks despite less stringent restrictions on the community, by using a refined strategy of targeted shutdown, restrictions on cases, their close contacts, identified hotspots and venues of concern rather than a whole of community lockdown. Risk assessments of potential transmission sites were constantly updated through our evolving experience with transmission events. However, this refined strategy did leave the potential for large point source outbreaks should any cases go undetected. [Addendum] A fourth wave that began in Sydney in June 2021 challenged this strategy due to the more transmissible nature of the Delta variant of SARS-CoV-2. LESSONS LEARNT: A wave of COVID-19 infections can develop quickly from one infected person. The community needs to remain vigilant, adhering to physical distancing measures, signing in to venues they visit, and getting tested if they have any symptoms. Signing out of venues on exit allows public health resources to be used more efficiently to respond to outbreaks.

An emergent clade of SARS-CoV-2 linked to returned travellers from Iran.

The SARS-CoV-2 epidemic has rapidly spread outside China with major outbreaks occurring in Italy, South Korea, and Iran. Phylogenetic analyses of whole-genome sequencing data identified a distinct SARS-CoV-2 clade linked to travellers returning from Iran to Australia and New Zealand. This study highlights potential viral diversity driving the epidemic in Iran, and underscores the power of rapid genome sequencing and public data sharing to improve the detection and management of emerging infectious diseases.

SARS-CoV-2 neutralizing antibodies: Longevity, breadth, and evasion by emerging viral variants.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) antibody neutralization response and its evasion by emerging viral variants and variant of concern (VOC) are unknown, but critical to understand reinfection risk and breakthrough infection following vaccination. Antibody immunoreactivity against SARS-CoV-2 antigens and Spike variants, inhibition of Spike-driven virus-cell fusion, and infectious SARS-CoV-2 neutralization were characterized in 807 serial samples from 233 reverse transcription polymerase chain reaction (RT-PCR)-confirmed Coronavirus Disease 2019 (COVID-19) individuals with detailed demographics and followed up to 7 months. A broad and sustained polyantigenic immunoreactivity against SARS-CoV-2 Spike, Membrane, and Nucleocapsid proteins, along with high viral neutralization, was associated with COVID-19 severity. A subgroup of "high responders" maintained high neutralizing responses over time, representing ideal convalescent plasma donors. Antibodies generated against SARS-CoV-2 during the first COVID-19 wave had reduced immunoreactivity and neutralization potency to emerging Spike variants and VOC. Accurate monitoring of SARS-CoV-2 antibody responses would be essential for selection of optimal responders and vaccine monitoring and design.

Risk factors leading to COVID-19 cases in a Sydney restaurant.

OBJECTIVE: To explore the factors associated with the transmission of SARS-CoV-2 to patrons of a restaurant. METHODS: A retrospective cohort design was undertaken, with spatial examination and genomic sequencing of cases. The cohort included all patrons who attended the restaurant on Saturday 25 July 2020. A case was identified as a person who tested positive to a validated specific Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) nucleic acid test. Associations were tested using chi-squared analysis of case versus non-case behaviours. RESULTS: Twenty cases were epidemiologically linked to exposure at the restaurant on 25 July 2020. All cases dined indoors. All cases able to be genomic sequenced were found to have the same unique mutational profile. Factors tested for an association to the outcome included attentiveness by staff, drink consumption, bathroom use and payment by credit card. No significant results were found. CONCLUSION: Indoor dining was identified as a key factor in SARS-CoV-2 transmission, and outdoor dining as a way to limit transmission. Implications for public health: This investigation provides empirical evidence to support public health policies regarding indoor dining.

Epidemiologic Evidence for Airborne Transmission of SARS-CoV-2 during Church Singing, Australia, 2020.

An outbreak of severe acute respiratory syndrome coronavirus 2 infection occurred among church attendees after an infectious chorister sang at multiple services. We detected 12 secondary case-patients. Video recordings of the services showed that case-patients were seated in the same section, up to 15 m from the primary case-patient, without close physical contact, suggesting airborne transmission.

Genetic diversity of SARS-CoV-2 and clinical, epidemiological characteristics of COVID-19 patients in Hanoi, Vietnam.

A second cluster of COVID-19 cases imported from Europe occured in Vietnam from early March 2020. We describe 44 SARS-CoV-2 RT-PCR positive patients (cycle threshold value <30) admitted to the National Hospital for Tropical Diseases in Hanoi between March 6 and April 15 2020. Whole SARS-CoV-2 genomes from these patients were sequenced using Illumina Miseq and analysed for common genetic variants and relationships to local and globally circulating strains. Results showed that 32 cases were Vietnamese with a median age of 37 years (range 15-74 years), and 23 were male. Most cases were acquired outside Vietnam, mainly from the UK (n = 15), other European countries (n = 14), Russia (n = 6) and countries in Asia (n = 3). No cases had travelled from China. Forty-one cases had symptoms at admission, typically dry cough (n = 36), fever (n = 20), sore throat (n = 14) and diarrhoea (n = 12). Hospitalisation was long with a median of 25 days, most commonly from 20-29 days. All SARS-CoV-2 genomes were similar (92-100% sequence homology) to the reference sequence Wuhan_1 (NC_045512), and 32 strains belonged to the B.1.1 lineage. The three most common variants were linked, and included C3037T, C14408T (nsp12: P323L) and A23403G (S: D614G) mutations. This group of mutations often accompanied variant C241T (39/44 genomes) or GGG 28881..28883 AAC (33/44 genomes). The prevalence of the former reflected probable European origin of viruses, and the transition D614G was dominant in Vietnam. New variants were identified; however, none could be associated with disease severity.